UTTG-18-93 



Effective Action and Renormalization Group Flow of 
Anisotropic Superconductors 



Steven Weinberg^ 

Theory Group, Department of Physics, University of Texas 
Austin, TX, 78712 



Abstract 

We calculate the effective action of a superconductor, without as- 
suming that either the electron-electron potential or the Fermi sur- 
face obey rotational invariance. This approach leads to the same gap 
equation and equilibrium free energy as more conventional methods. 
The results are used to obtain the Gell-Mann - Low renormalization 
group equations for the electron-electron potential. 
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1. Introduction 

This paper aims to demonstrate the usefulness in the theory of supercon- 
ductivity of the effective action formalism of quantum field theory^. Although 
the effective action may be defined non-perturbatively, for our purposes it 
may be taken simply as the sum of all connected one-particle-irreducible 
vacuum diagrams in the presence of any background field. As we shall see, 
the usual assumptions of the BCS superconductivity theory allow an almost 
trivial calculation of the effective action, from which one may obtain all of 
the properties of the superconductor, including the gap field, free energy, 
penetration depth, and so on. 

As an additional application of this formalism, we shall derive the renor- 
malization group flow of the electron-electron interaction, a subject that has 
received increased attention in the last few years [2]. The renormalization 
group equation for the electron-electron potential is obtained here from the 
condition that the effective action expressed as a function of the gap field 
should be renormalization group invariant. From a practical point of view, 
the most significant difference with most earlier work is that here we make 
no special assumptions about the form of the Fermi surface or the rotational 

invariance of the electron-electron potential.^ Just as in the rotationally in- 

2 For a lucid review of the effective action method and references to earlier work, see S. 
Coleman, ref. [1]. 

3 After this work was completed I learned that the renormalization group equation for 
non-circular Fermi surfaces in two dimensions had been derived by R. Shankar, Rev. Mod. 
Phys., ref. [2]. 
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variant case, it turns out that the renormalization group equations may be 
expressed as de-coupled equations for an infinite number of coupling param- 
eters. 

2. Calculation of the Effective Action 

As usual in the Fermi theory of liquids, we assume that all degrees of free- 
dom may be "integrated out" except for electrons (strictly speaking, quasi- 
particles with the quantum numbers of electrons) in a narrow shell around a 
Fermi surface, say of thickness \i in electron energy. (In the approach to be 
used here, this condition on electron wave numbers is implemented through 
the renormalization procedure, and must be verified a posteriori.) Also as 
usual, we discard all interactions that become irrelevant in the limit /i — > 0. 
This leaves just the one-electron and two-electron terms in the Lagrangian: 



where is an external electromagnetic vector potential, and s, s', etc. are 



adjusted to be a simple Kronecker delta by suitable definition of the electron 
field operator ip(x,s). The time argument is suppressed in all field opera- 
tors.] Also, as well known, the only diagrams that are not suppressed by 
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i— + eA (x, t) 5 s i tS + E s > >s (-iV + eA(x, t)) ip(x, s) 



x ^(x[, s[) ij)\xf 2 , s' 2 ) ip(x 1: Si) ip(x 2 , s 2 ) 



(1) 



spin indices summed over the values ±|. 



[The coefficient of the first term is 



powers of \i as [i — > are those whose structure constrains each interacting 
pair of electrons to have opposite momenta (for very slowly varying exter- 
nal fields), so that if one electron is on the Fermi surface, then the other is 
also. [The Fermi surface is defined by the vanishing of any eigenvalue of the 
energy matrix E s / jS (p), which is understood to include a chemical potential 
term. Time-reversal invariance tells us that if a momentum p is on the Fermi 
surface, then so is —p.] In particular, the unsuppressed vacuum diagrams 
are those that become disconnected if we cut through an interaction vertex 
so as to separate the incoming from the outgoing electrons. 

Before setting out to calculate these diagrams, we will introduce a pair 
field \I/ by a familiar trick, generally attributed to Hubbard and Stratonovich[3], 
which we extend here to general potentials. Add a term to the Lagrangian 
of the form 



and integrate over the pair field \l/(aTi,a^) as well as over the electron field 
ip(x). This clearly has no effect; the action is quadratic in \1>, so the path 
integral may be evaluated by setting \1> equal to its equilibrium value 

^ sl)S2 (xi,£ 2 ) = ^(fi, Si) lp(x 2 , S 2 ) 

at which (2) vanishes. Instead of integrating over the pair field, we shall 




x Oci,^) - i/>(x u si) ip(x 2 , s 2 )} 



(2) 
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evaluate the effective action in the presence of a background pair field, inte- 
grating over the electron field. The term (2) has been chosen to cancel the 
term in the Lagrangian quartic in the electron field, leaving just quadratic 
terms: 



L + AL = Y / Jd 3 x ^(x, s') (-i^ + eAo(x, t)^j 5 S ,, S + E s , yS (-iV + eA(x, t)) ip(x, s) 
-- J d 3 x[d 3 x' 2 d 3 x 1 d 3 x 2 V s >^ SUS2 {x\\x 2 \x\,x 2 ) 

s' 1 ,s' 2 ,si,S2 

x [^(^i, si) ^( x 2, 4) ^LoaO^ij ^) + ^(^1, «i) ^{x 2 , s 2 ) (a^, x' 2 ) 

+ 4 S J d 3 x' 1 d 3 x' 2 d 3 x 1 d 3 x 2 ^^(xi, 4) K i)S ^ Slifl2 (xl',x^',xl,^) # sl , S2 (aTi, aT 2 ) 

(3) 



S, ,So,Sl,S2 



Now, the effective action in a background pair field is given by the sum of 
all one-particle-irreducible vacuum diagrams — that is, all vacuum diagrams 
that cannot be disconnected by cutting through any one internal line[l]. On 
the other hand, we have already mentioned that in the limit \i — > 0, we are to 
keep only graphs that can be disconnected by slicing through any electron- 
electron interaction so as to separate incoming from outgoing electrons. In 
using (3), this means that we are to keep only graphs that are disconnected by 
cutting through any internal pair field line. Since we are to keep only graphs 
that both are and are not disconnected by cutting through any internal pair 
field line, we conclude that we must keep only graphs that have no internal 
pair field lines at all. There are just two such graphs; a tree graph arising 



from the last term in (3), and a one-electron-loop graph whose value is given 
by the determinant of the "matrix" accompanying the terms in (3) quadratic 
in the electron field: 

r[tf] = - I dt I d 3 x[d 3 x' 2 d 3 x 1 d 3 x 2 



s', ,s',s ly s 2 



— In Det 
2 



A B 

-A T 



+ - In Det 
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A 
-A T 



+ r[o] (4) 



where A and B are the "matrices" : 



A - = 



d \ - 

-i — + eA (x, t) J 6,,,, + E s , jS (-iV + eA(x, t)) 



5 3 (x' - x)5(t' 



®s> X -t>, sS t = ^s's{x>,x)5{t'-t) 

and A s > s (x',x) is the "gap" field: 

A s , s (x>, x) = ~Y,j d3 y d3 y' v s>s*'<j(x', x, y', y)^ a , a {y'y) (7) 

<r' ,a 

The constant T(0) represents the contribution of electrons not near the Fermi 
surface, for which the approximations made here are not applicable. 

To avoid becoming lost in a cloud of indices, we now specialize to the case 
of spin-independent forces and a spin-singlet pair field, writing 



vj> 1 _i 

2 2 



1 = * 

2 2 



Vj/l 1 

2 2 



2 2 



Vi _i i _i — Vi _i _i i 

2222 2222 



— V_l 1 1 _1 + V_l 1_ 1 I 

2222 2222 



E5 s t s 



(8) 

2V (9) 
(10) 



(It is easy to extend the results here to systems with a more general spin 
dependence, such as liquid He 3 .) It follows from (8) and (9) that 



Ai _i = -A_i i = A An = A_i _i = 

2 2 22 22 22 

where 

A(x', x) = - f d 3 yd 3 y' V(x', x, y f , y) ^{y'y) 
The effective action is now 



(11) 



(12) 



r[tf] = J dtf d 3 x' 1 d 3 x' 2 d 3 x 1 d 3 x 2 V{xi,X2,x u x 2 ) *(x*i,x 2 ) 

A 



—i In Det 



+ i In Det 



A B 

& -A T ^ " ^' -A 1 
d 

-i— + eA {x, t) + E{-iV + eA(x, t)) 



+ r[o] 



(13) 



5 3 (x'-x)5(t' -t) 



Bp-» = A(:?,iO*(f-f) 



(14) 
(15) 



3. Translationally Invariant Equilibrium 

First we consider the translationally invariant case, with no external elec- 
tromagnetic fields. Then the pair and gap fields can be put in the form 

ty(x>,x) = J d^pe^'-^ip) (16) 
A(x', x) = (2n)' 3 f d 3 p A(p) (17) 



The electron-electron potential appears here in the form 

J d 3 Xl d 3 x 2 d 3 x 3 d 3 x 4 e^*-*"*) e ^(*s-i4) V{x 1 ,x 2 , x 3 , x 4 ) 
= Q,V(p',p) 



(18) 



where f2 4 is the spacetime volume 

tt 4 = I d 3 x I dt (19) 



The effective potential U[$] is definedfl] as minus the effective action per 
spacetime volume 

um = -r[*]/n 4 

= -J d 3 pd 3 p' y*(p')V(p',p)^(p) 

|A(p)p 



^/^H 1 - ^^ )'^ 01 (20) 



with 



A(p) = - J d 3 p V(p,p') (21) 

(We are working here at zero temperature. For non-zero temperature, the 
integral over uj is of course replaced with a sum over discrete Matsubara 
frequencies.) Wick rotating, integrating over cu, and expressing \1> in terms 
of A, this becomes 

U[A] = - J d 3 pd 3 p' A*(p')V-\p',p)A(p) 

We pause to note that the gap equation is obtained from the condition 
that A(p) be at the stationary point A (p) of U[A]: 
5U[A] 

A=A 



= 



SA*(p) 



= -[ d 3 p' V-\p,p')A (p') - -j-- . AoW 

J 2(2^)3/^(^ + 1^^12 



or in a more familiar form 



J3 K(p,j/)A (p') 



2 ( 27r ) 7 ^(pO + IAo^OI 2 



Also, the effective potential at this stationary "point" is the free energy 
density, which using (23) may be expressed in terms of the gap field: 



^equilibrium ~ Fa=0 — U[Ao] — U[0] 

lAobll 2 



(2tt)* 



d 6 p 



E\p) + \A Q (p)\i + E(p) 



[24) 



_2^EP(0) + \A o (p)\* 
This is the same as the result normally derived from the BCS ground state 
wave function [4]. However, it should be noted that (22) is not the same as 
the formula for the free energy as a functional of A(p) calculated from the 
BCS wave function^: 

= +J d 3 pd 3 p' V(p',p) *(p) + J d 3 pE{p) 



1 - yjl -4(2tt) 6 |^ 
(25) 

This difference arises because the effective potential does not have the inter- 
pretation of an energy density for a composite field like *&(p) or A(p), except 
at the stationary point of the effective potential[6]. Nevertheless, both (22) 
and (25) yield the same gap equation, and the same value for the equilibrium 
free energy. 

4 This may be obtained e. g. from Eq. (4-64) of ref. [5], using the relations = ut Vk 
and u| + v\ = 1 to express the coefficients Uk and in terms of the pair field , and 
then converting to the notation of the present paper by replacing sums over the discrete 
index k with integrals over p, and inserting appropriate factors of (27r) 3 . 
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So far, the limitation of electron momenta to a thin shell around the 
Fermi surface has been left implicit. To make this limitation explicit, we will 
write the electron momenta as 

p = k + h(k)l d 3 p = d 2 kd£ (26) 

where k is on the Fermi surface (that is, E{k) = 0), and h{k) is the unit vector 
normal to the Fermi surface at k. [Here d 2 k should be understood as Jd9 1 d6 2 , 
where Q\ and 6 2 are coordinates on the Fermi surface, and J is the Jacobian 
of the transformation from p to £, 6>i, and 6> 2 . For a spherical Fermi surface, 
/ d 2 k just gives a factor 4irk F .] We will temporarily impose the condition 
that electron momenta are close to the Fermi surface by introducing a cut- 
off A on £, chosen small enough so that A <C \k\ for all k on the Fermi surface, 
and so that V(p',p) and A(p) change negligibly as i and i' vary from to A. 
[This cut-off will eventually be obviated by the introduction of a renormalized 
electron-electron potential.] With this cut-off, we may approximate 

E(p) = v F (k)£ (27) 

where v F is the Fermi velocity: 

(28) 



, A dE{k + n(k)£) 
v F (k) = — 



dl 

Eq. (22) may now be written 

U[A] = -A 2 J d 2 kd 2 k' A*(k')V~ 1 (k', k)A(k) 



1 r A 

' di I d l k 



(2vr) 3 Jo 



£ 2 v 2 F (k) + \A(k)\ 2 - £v F (k) + U[0] (29) 



so that U[A] is now denned as a functional of the gap field on the Fermi 
surface only. We define a renormalized electron-electron potential at a renor- 
malization scale \x by 

5 2 U[A] 



V~\hf,k) = - 



(30) 



A(Jfc)=A(fc)*=/i 



5A*{k')5A(k) 

When we express the electron-electron potential V in terms of the renormal 
ized potential V^, Eq. (29) for the effective potential becomes: 

U[A] = - J d 2 kd 2 k' A*(k>)V~\k>,k)A(k) 

^j Q A d£ J d 2 k \jl 2 v 2 F {k) + |A(fc)| 2 - £v F (k) 

„ 2 \A(h\\ 2 

+ U[0] 



|A(*)| ! 



+ 



(31) 



2{Pv 2 F {k) + // 2 )V2 4(£ 2 Vp(k) + /i 2 ) 3 / 2 
The one- loop integral over i now converges as I — > oo, so we may remove the 
cut-off, and find 



[/[A] = ~y c^Jfed 2 *;' A\k')V-\k',k)A{k) 



2(2tt) ; 



|A(AQ| a 



In 



- 1 



+ u[o] 



(32) 



The corresponding gap equation is obtained from the condition that this 
expression for the effective potential be stationary: 



Ao(fc) 



1 



2(2tt) 3 

and the equilibrium free energy is 



d 2 k'V fl (k,k') v F \k')A {k') 



In 



|A (fcQ| 



(33) 



equilibrium 



Fa=o — 



1 



4(2tt) 3 
10 



d 2 k 



|A (fc)p 
v F (k) 



(34) 



These very simple results should not be taken entirely literally, because 
the effective potential will always contain "irrelevant" terms of higher order 
in A arising from degrees of freedom that have been integrated out here. In 
particular, note that according to (32) the difference U[A] — U[0] vanishes 
for A — > 0, is negative for sufficiently small A, and goes to +00 for A — > 00, 
so it must have a stationary point A 7^ 0, at which the gap equation (33) is 
satisfied, for any renormalized electron-electron potential, whether repulsive 
or attractive. But this solution should not be taken seriously if it has A so 
large that it is outside the range of validity of these equations. To be specific, 
it is easy to see from (32) that if V^(k',k) is positive (in the matrix sense) 
for some /z, then the minimum A (A;) of the effective potential will have a 
scale ||A || > y/e/j,, where the "scale" ||A|| of an arbitrary function A(/c) is 
defined by the condition 



frequency, then the solution of the gap equation (33) is physically irrelevant. 

4. Renormalization Group Flow 

In the Wilson approach to the renormalization group that is most familiar 
in condensed matter physics, we would keep the cut-off A finite, and derive a 
renormalization group equation for the A-dependence of the unrenormalized 
electron-electron potential V(k', k) from the condition that the effective po- 




In particular, if V^(k', k) is a positive "matrix" for fi of the order of the Debye 
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tential (29) should be A-independent. This approach would be possible here, 
but it would require the introduction of "irrelevant" terms in U[A] of higher 
order in A(k) to keep U[A] A-independent for finite A. It is much simpler to 
apply the older approach of Gell-Mann and Low: the renormalization group 
equation in this approach is the condition that U[A] is independent of the 
arbitrary renormalization scale \i: 

^ V ^'^ = - Sf'~% (35) 
dp " 2(2ir) 3 v F (k) 

or equivalently 

* T/^'' l) = 2(ky I d2k " v »$' g,) v ^ 1{kl/) v;(A;7, ' k) (36) 

This of course also implies that the solution of the gap equation is inde- 
pendent of fx. For the special case of a spherical Fermi surface Eq. (36) is a 
continuous version of the discrete renormalization group equation of Benfatto 
and Gallavotti[2], provided we identify their constant (3 with 1/2(27t)Vf- Eq. 
(36) also agrees with the results of Shankar[2], with fairly obvious changes 
to convert his results from two to three dimensions. 

These renormalization group equations can be separated into equations 
for the eigenvalues A„(/i) of the Hermitian kernel 

K„(k>, k) ee v- 1/2 {k>) v~ 1/2 (k) V,(k>, k) . (37) 

From either (35) or (36), we see that the eigenvectors of K are renormalizat ion- 
group invariant, while the eigenvalues are governed by the flow equations 

^n(^) = KM ■ (38) 
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For a spherical Fermi surface these eigenvectors are just the spherical har- 
monics Yg n (k), but we see that the decoupling of the eigenmodes of K is in 
fact quite general, not depending on rotational invariance. 

A completely repulsive potential may be defined as one for which all eigen- 
vectors A n are positive. If this is true at some starting scale /j,q (say, the Debye 
frequency) then as ji — > the eigenvalues stay positive and become smaller, so 
nothing interesting happens. This conclusion may be altered by the inclusion 
of formally irrelevant electron-electron couplings, as in the Kohn-Luttinger 
effect, but is not directly affected by anisotropies. A completely or partly 
attractive potential is one for which all or some of the eigenvalues of K are 
negative. Any eigenvalue that is negative at some /x becomes infinite at 
a finite ji < /x , with the eigenvalue that is largest in magnitude becoming 
infinite first. As already mentioned, this is the case where superconductivity 
actually occurs. 

It is not clear what is gained by this renormalization group analysis. 
Within Fermi liquid theory, the effective potential arises solely from tree 
and one-loop graphs, so we do not need to use the renormalization group 
to make the sort of improvement in perturbation theory familiar in quan- 
tum electrodynamics, or to identify a weak coupling regime, as in quantum 
chromodynamics or the theory of critical phenomena. Of course, one can go 
beyond tree and one-loop graphs, but this would require that we take into 
account irrelevant terms in the original Lagrangian, which would involve ad- 
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ditional microscopic information, not just renormalization group equations. 

5. Slowly Varying Electromagnetic and Goldstone Fields 

Now let us return to the general case of a superconductor in a translationally- 
non-invariant external electromagnetic field. In the limit where this field has 
very small frequencies and wave numbers (smaller than the inverse correla- 
tion length), we can integrate out all degrees of freedom except those that 
have zero "mass," in the sense that their frequency vanishes in the limit of 
vanishing wave number. For a superconductor that is not close to a phase 
transition at which superconductivity is lost, the only such "massless" degree 
of freedom is the Goldstone mode, associated with the spontaneous break- 
down of electromagnetic gauge invariance within the superconductor. All of 
the classic exact properties of superconductors (persistent currents, Meissner 
effect, flux quantization, and Josephson frequency) can be derived by consid- 
ering only general properties of the effective action for the Goldstone mode 
in the presence of external electromagnetic fields [7]. But to derive values for 
quantities like the penetration depth in Type II superconductors, we need a 
detailed formula for the effective action. 

The effective action for a Goldstone mode is in general obtained by setting 
all the fields in the effective action equal to their equilibrium values, and then 
subjecting them to a symmetry transformation with space-time dependent 
parameters equal to the Goldstone fields. In our case, where the broken 



14 



symmetry is electromagnetic gauge invariance, this means that we must make 
the replacement 



A{x', x, t) -> e-^^Aoix' - x)e~ ie ^ 



(39) 



where (f>(x) is the Goldstone field, and A (x' — x) is the equilibrium gap field, 
given by the Fourier transform of the solution of the gap equation: 

A (x> -x) = (2tt)- 3 J d 3 p e^'-^Aoip) (40) 

Electromagnetic gauge invar iancef] then allows us to remove the factors e~ ie< t>( x ' ) 
and e~ ie ^ x '^ in Eq. (39), by replacing the electromagnetic vector potential 
Afj,(x, t) with Afj,(x, t) — d^(p(x, t). In this way, Eqs. (13)-(15) yield the effec- 
tive action: 

'A 



V[<j>,A] = T A=0 [A] -i In Det 
where now 

d 



A B 

B+ -A T 



+ % In Det 



-A 2 



(41) 



-i— + eA (x, t) - ecj)(x, t) + E(-iV + eA(x, t) - eV0(f , t)) 



A- 

x 5 3 (x' -x)6(t' -t) 
B xM3t = A (x>-x)6(t'-t) 



(42) 
(43) 



Quantitative properties of the superconductor such as the penetration depth 
can be read off from the expansion of Eq. (41) in powers of A (x, t) — <f)(x, t) 

and A(x, t) — V0(x, t). 

5 The Lagrangian (1) is gauge invariant for either a local electron-electron potential, 
or an arbitrary electron-electron potential with a suitable dependence on the external 
electromagnetic fields. 
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